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Artificial biological pacemakers were developed and tested in canine ventricles. Next steps will 
require obtaining oscillations sensitive to external regulations, and robust with respect to long term 
drifts of expression levels of pacemaker currents and gap junctions. We introduce mathematical 
models intended to be used in parallel with the experiments. 

The models describe human mesenchymal stem cells {hMSC ) transfected with HCN2 genes and 
connected to myocytes. They are intended to mimic experiments with oscillation induction in a cell 
pair, in cell culture and in the cardiac tissue. We give examples of oscillations in a cell pair, in a one 
dimensional cell culture and dependence of oscillation on number of pacemaker channels per cell and 
number of gap junctions. The models permit to mimic experiments with levels of gene expressions 
that are not achieved yet and to predict if the work to achieve this levels will significantly increase 
the quality of oscillations. This give arguments for selecting the directions of the experimental work. 



PACS numbers: 

Introduction 

Genetically engineered pacemakers could be developed 
as a possible alternative to implantable electronic de- 
vices. Several approaches are developing: heterologous 
h2-adrenergic receptor expression in pig atria over- 
expression of If Q and dominant negative suppression 
of Iki in guinea pig ventricle Q. 

Impressive results were obtained with stem cells to cre- 
ate an artificial biological pacemaker Q: human mes- 
enchymal stem cells {hMSC) transfected with HCN2 
genes injected in canine left bundle branch provide ven- 
tricular escape rhythms that were mapped to the site 
of the injection. This experimental success shows great 
promises for future therapy. 

In a recent review article 5], the question : "What 
characteristics should be embodied in a biological pace- 
maker?" was put forward. The creation of a stable phys- 
iological rhythm was listed as one of the most important 
properties for future biopacemakers. The present article 
describes the mathematical models intended to be runned 
in parallel with the experiments to achieve this goal. 

The oscillations should be reasonably stable with re- 
spect to inevitable variations of the parameters affecting 
the biological pacemaker. Indeed, many factors can af- 
fect the period of the pacemaker namely: decreasing the 
expression levels of pacemaker channels or gap junctions, 
modification of the peptides regulating gene expressions, 
loss of genes, migration of small percentage of stem cells 
to other sites or differentiation into other cell types; loss 
of pacemaker function by some cells. 

In order to obtain a stable pacemaker, many time con- 
suming experiments should be performed. In exploring 
such systems, we believe that an adequate combination 
of experiments with theoretical predictions can lead to 
a much faster results, by reaching a better quantitative 
understanding and in so doing, in reducing the number of 



experiments. To this end, one needs precise and adapt- 
able mathematical models based on well-established con- 
cepts and results |6|-p3|. 

The models presented here describe human mesenchy- 
mal stem cells {hMSC ) transfected with HCN2 genes and 
connected to myocytes. They are supposed to be runned 
on a computer mimicking the planned experiments. A 
part of these computer simulations can be used instead 
of the planned experiments while experiments performed 
in parallel with the computer simulations can be used to 
update the models. 

In the first version of the models, all the complex pro- 
cesses mentioned above are condensed into three control 
parameters: the number of pacemaker channels per stem 
cell, N] the number of gap junctions per cell, n and the 
concentration of the stem cells in the pacemaker. 

We demonstrate how they work by mimicking exper- 
iments with oscillations induction in a cell pair, in cell 
culture and in the cardiac tissue intended to study the 
robustness of oscillations and stability of the oscillation 
period. 

The models permit to easily mimic experiments with 
arbitrary levels of gene expressions (including levels not 
yet achieved in experiments) and to assess whether 
achieving such levels of expression will increase the qual- 
ity of oscillations. 

Models 

Stem cell models 

Models for a stem cell hMSC transfected with HCN2 
gene were constructed from experimental data Q • We 
focus here on two models: a model, including only in- 
activation similar to models of HCN channels by [2ll| . 
0,23 1 and a model with two gating variables, acti- 
vation and inactivation of the pacemaker current, similar 
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to the Hodgkin-Huxley (HH) model [3. 
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FIG. 1; Stem cell model, a) steady values of activation g\ 
and inactivation g-2,. b) time constants. 

Stem cell model with inactivation only (Stem-h) is: 

dEs/dt - -Naf 5| [Es - Er]/Cs (la) 
d9'2/dt^[r2{Es)-9'2]/rm) (lb) 

Here, Eg is the stem cell membrane potential, mV. N 
is the number of pacemaker channels per cell, tx/ ~ IpS 
is the conductivity of a single pacemaker channel, C'^95 
pF is the stem cell capacitance, iJ^ ~ - 40 mV is the 
reversal potential of the HCN2 pacemaker current. 

g| is a gating variable describing inactivation kinetics 
(notation h is used in Hodgkin-Huxley eqns'Ts'l). 5|(i?s) 
is its steady value that depends on Eg, r|(i?s) is the 
inactivation time constant that depends on Eg as well. 

The model (Stem-h) is totally defined by two functions 
only, ^|(i?s), T2{Es); their graphs are shown in Fig^ 

Comment 1. It is possible to use more detailed nota- 
tions keeping track that the pacemaking current It con- 
sists of two components, I^a and Ik (e.g., as in [l3l[l5| ). 

dEg/dt = If^Na + Ifji (2a) 

If,Na = gf.Na Y [Eg - EMa] (2b) 

If.K = gf.K Y [Eg ~ Ek] (2c) 
dY/dt = {Y'{Eg) - Y)/T^{Eg) (2d) 



But since both components of the current pass through 
the same ionic channel and thus are gated by the same 
gating variable Y, this description collapses back to Eq^ 
with 

Er — [En a 9f,Na + Ek 9f,K]/[gf,K + gf,Na] (3) 

Comment 2. A model with one gating variable only 
(slow inactivation) is used in all descriptions of the HCN 
currents that we know. The justification of the above 
statement is that the pacemaker current induces oscilla- 
tions with period 1 sec and thus there is no need for 
description of processes on a faster time scale. 

But this is not a reliable argument. The front of an 
AP induced by HCN develops on a much shorter time 
scale. To accurately describe HCN effects on it, the 
much shorter processes should be taken into account. 
So, we developed also a model incorporating two gating 
variables, fast activation and slow inactivation (similar 
to m and h in HH equations). 

Stem cell model with activation and inactivation 
(Stem-mh) is: 

dEg/dt = ^Nafglg^Eg - Er)/Cg (4a) 
dg!/dt^igt-gt)MiEg), z = l,2 (4b) 

Eg is the stem cell membrane potential, mV. gf, g| are 
activation and inactivation gating variables (notation 
m,h are used in Hodgkin-Huxley eonsfl^'). T;f,T| are 
activation and inactivation time constants, Fig^ 

Different types of myocytes 

We need to investigate interaction of stem cells with 
myocytes of different types. No published models are 
sufficient for this (every model describes a specific type of 
myocyte) . For this reason, we will content ourselves with 
an existing model, which we will modify to reproduce 
qualitatively the various types of myocytes. 

In this article, we work with several types of myocyte: 
ventricular adult myocyte (with high threshold logc for 
inducing oscillations); self-oscillating neonatal type my- 
ocyte and several transitional types. We obtain mod- 
els for these myocytes from a model for ventricular my- 
ocytes, using different values of expression levels of Iki 
current and re-scaling the time constants to keep the ac- 
tion potential duration (APD) identical for all myocytes. 

So, Iki is expressed as: 

Iki = k' l},, (5) 

where k' is the number of Iki channels per myocyte 
and I^^ is the current through one Iki channel (averaged 
in time). 

Our models were derived from the Beeler-Reuter model 
(BR) which describes the membrane potential, Em, 
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gating variables controlling ionic chanels and Calcium 
concentration [Ca\: 

dtEm = -hon{Ern,gr)/Crn+DV^E„, (6a) 

dtgr = {-g? - gT)lTT{E^). * = 1,6 (6b) 

dt\Ca\ = -lO^^ica + 0.07(10"'^ - [Ca]) (6c) 

Here ica = gcag^ge^i^rn- Eca) , gca = OMmS/cm^ 
and Eca = -82.3 - 13.0287 ln[Ca]. The diffusion term 
D in EqEl represents the electric coupling between cells: 
D ~ lQ~^c'm? /raSec. The model was calculated in one 
and two dimensional cases with dx = O.lmmanddt — 
OmmSec. 

A myocyte with k'— 1400 Iki channels per myocyte 
corresponds to the original BR model. We will use 
also dimensionless values k = fc'/1400. The excitation 
threshold is lexc = 2pA/pF and the oscillation threshold 
lose = 2.01pA/pF. The myocytes are auto-oscillatory for 
k' <~ 320 (fc <~ 0.23) and are excitable for k' 320. 

When k' approaches from below to the critical value 
fc^ ~ 320, the amplitude of the oscillation remains the 
same, while the period T^sc of the oscillation diverges 
logarithmically Tosc ~ ln{kc — k), k < kc- This indicates 
the change from oscillatory to excitable but quiescent be- 
havior occurs through a homoclinic bifurcation [23|.[T^'). 
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FIG. 2: Decrease of oscillations period in neonatal 
myocytes by HCN2 transfected stem cells. A cell cul- 
ture consisting of neonatal self-oscillating cells was simulated. 
A record far from the pacemaker region: M- myocytes alone; 
MS- stem cells are added. A one dimensional simulation: 
201 myocytes with k'— 300 Iki channels per myocyte, my- 
ocytes alone oscillate with period Tb=650ms. Stem cells with 
N = 10^ pacemaker channels per stem cell, n= 1000 gap 
junctions are added to the pacemaker region I = 4.3mm long. 
They decreased the period to 425 ms. 

Cell pair model 

A " cell pair" , consisting of a myocyte and a stem cell 
coupled together is described by the following set of equa- 
tions : 
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FIG. 3: Induction of oscillations by a stem cell con- 
nected to a myocyte. Myocyte and stem cell potentials 
are noted as M and S. Above each image the number of gap 
junctions n per cell is shown, a) cells are not connected (n=0), 
the myocyte generates one AP only, b) oscillations in the cell 
pair, c) disappearance of oscillations when the cell connection 
is too strong. Parameters: the number of pacemaker channels 
per stem cell A'^ = 1.2 x 10^ (see eq|7^), the number of Iki 
channels per myocyte k' — 326 (see eql^. 



dtEs = -Nafg'2{Es+iO)/Cs + 

nag{Em ~ Es)/Cs (7a) 

dtEm — —Iion(Em, gT") /Cm -\- 

nGg{Es-Era)/Cra (7b) 

dtg? = F{gr,Em)/Tr{Em), i = l,...,6 (7c) 
dt[Ca] ^ -lO^^ica + 0.07(10"^^ - [Ca]) (7d) 
dg'Jdt = F{gi,E,)/T^{E,) (7e) 

where n is the number of gap junction channels per cell, 
Ug is the conductivity of a single gap junction channel 
(50 pS). Term nag{Em — Eg) in Eq[7^ describes current 
flowing from a myocyte to a stem cell. 

Here, we described the stem cell by the simplified 
model EqQ The implementation of the full model Eq01 
leads to a very similar system. 

Cardiac tissue with incorporated stem cells 

We suppose that stem cells are injected to a pacemaker 
region only and are connected to myocytes. The tissue 
is described by EqsElEI Inside the pacemaker region, to 
Eq|S^, the term nag{Es — £",„) is added, similar to Eq. 
{7}p. Numerically, we considered both a one dimensional 
(ID) fiber and a two dimensional (2D) tissue. 

Increase of oscillation frequency 
in self-oscillatory myocytes. 
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Gap junction channels/Cell 

FIG. 4: Oscillations phase plane for a cell pair. Bound- 
aries of oscillation regions for two models are shown: a model 
with inactivation only (model Stem-h, Eq^ and a model with 
activation and inactivation (model Stem-mh, Eq|lJ. Note 
that in both models, high expression levels of gap junctions 
or pacemaker channels are not needed to induce oscillations 
and may even destroy them. 



In experiments '6| , stem cells added to cultured neona- 
tal myocytes that are self-oscillatory. The stem cells 
decreased the oscillation period of these myocytes from 
T ^ 660 to r ^ 370ms. This phenomenon is easily repro- 
duced by the model, both in a cell pair and in a cardiac 
tissue, FigEl 

This decreasing of oscillation frequency by HCN genes 
is a robust phenomenon. In addition, it can be observed 
at arbitrary values of parameters; only the amount of the 
period decrease is parameter dependent. 

A large interval of number of pacemaker channels per 
cell TV and number of gap junction channels per cell n 
can be easily investigated. In todays experiments, only 
a part of this interval is accessible. 

Induction of oscillations in quies- 
cent myocytes. 

Induction of oscillations in a cell pair 

Induction of oscillations by a stem cell connected to 
a myocyte is shown in FigEl When the cells are not 
connected (Fig|3t, n=0), the myocyte can generate one 
AP only. When the cells are connected (Figl^D, n=20), 
a pacemaker activity can be initiated in the cell pair. A 
diastolic depolarization is observed in the myocytes, and 
it generates large amplitude periodic APs. 

This result confirms the conjecture, based on results of 
experiments with tissues 0, [qI, that a stem cell hMSC 
transfected with HCN2 gene can induce oscillations in a 
cardiac cell. 

Disappearance of oscillations 

When the cells are too strongly connected (FigEt, 
n=200), the pacemaker activity decays: the diastolic po- 



tential shifts closer to the resting potential of the stem 
cell {Er ^ — 40mT^) and the sodium channels are inacti- 
vated. This is not a sharp termination of oscillations: the 
oscillation amplitude decreases smoothly as the number 
of gap junctions n increases. When the amplitude of gen- 
erated periodical AP's is too small, they will be unable 
to propagate in cardiac tissue. 

In Fig01 the position of " disappearance of oscillation" 
boundary was chosen at parameter values where the up- 
stroke of the AP decreases below ~ mV and the ampli- 
tude of oscillations becomes inferior ^ 50 mV. It would 
be more accurate to determine when the propagation will 
disappear using cardiac tissue models. 

Robustness of oscillations. 

As we observed from the example above, oscillations 
in excitable cardiac cells can be induced not for all pa- 
rameter values. 

In Fig^lphase plane for a cell oscillation pair is shown. 
It is seen that oscillations are not observed at low values 
of number N of pacemaker channels per cell and num- 
ber n of gap junction channels per cell. When n, N are 
increasing, oscillations arrive as a jump, with full ampli- 
tude but with a very high oscillations period (several sec- 
onds). The period decreases and becomes normal when 
N,n continue increasing. It is the same homoclinic bifur- 
cation 24] that describes the transition from oscillatory 
to excitable but quiescent behavior when expression level 
of Iki is increased. 

At the contrary, when oscillations disappear in a cell 
pair, there is no jump: amplitude of oscillations slowly 
decrease while values of n, N increasing too large. 

Fig01 shows the oscillation boundaries for two models: 
the model comprising inactivation gating variable only as 
in majority of the descriptions published (model Stem-h), 
and the model comprising both gating variables, activa- 
tion and inactivation (model Stem-mh). Since the ratio 
e — T;f/r| ^ 1 (Fig. ^), it is natural to expect that a 
simplified model with inactivation only will give results 
close to the full model. As illustrated in Fig^the results 
are qualitatively but not quantitatively similar. 

As seen in Fig0] the induced oscillations are robust: 
they exist for a wide range of parameters: number of gap 
junctions per cell n and number of pacemaker channels 
per cell N. This is important to consider for constructing 
a biological pacemaker. 

Cell pair experiments. Induction of oscillations in a 
cell pair experiment was not demonstrated yet. What is 
the reason? A possible explanation : Expression of gap 
junctions in a cell pair 100 gap junctions/cell) is 1-2 
orders of magnitude less than in a normal cardiac tissue. 
Then, a path to induce oscillations in a cell pair might 
be increasing the expression of gap junctions. 

As Fig^ show, approximately 100 gap junctions per 
cell is well enough for induction of oscillations. The mod- 
els show that : 

- an experimental path intended to increase the expres- 
sion of gap junctions will not give positive results. 
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- selecting myocytes with not high level of expression 
of Iki will result in oscillations. 
Conclusions 

1. Mathematical models of stem cells (hMSC) trans- 
fected with HCN2 genes and connected to myocytes were 
developed. This opens a way to investigate the induction 
of oscillation and its robustness in a cell pair, in cell cul- 
ture and in the three dimensional cardiac tissue. 

2. The mathematical results confirm the conjecture, 
based on experiments with tissues, that stem cells can 
induce oscillation even in myocytes that do not self- os- 
cillate. The oscillations are robust and well suited for 
creation of biological pacemakers (Fig0J. 



3. To induce oscillations in a cell pair, there is no need 
to increase expression levels of gap junctions or of HCN 
genes (Fig0J. 

4. The models have a potential to eliminate a part 
of time consuming genetical and electrophysiological ex- 
periments needed for developing a biological pacemaker. 
They also permit to estimate the effects of expression 
levels that are not achieved in the experiments yet (Figs, 
ij. 
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